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Abstract. We explicitly solve the Goryachev case of rigid body motion in terms of sigma- functions of a 
trigonal curve of genus 3. Although completely integrable this system is not algebraically integrable, and 
the solution involves dynamics on the strata of the curve. The solution is not meromorphic in complex time. 
Although the strata of hyperelliptic curves have appeared in the literature in the context of some dynamical 
systems, this is the first example of an integrable system whose solution involves a more general curve and 
several new features (and formulae) are encountered in the solution. 



1. Introduction 



Most of the known finite-dimensional integrable systems of classical mechanics and mathematical physics 
are also algebraically completely integrable: following [AvM989] . their invariant tori can be extended to 
specific complex tori, Abelian varieties, and the complexified flow is a straight-line flow on them. As a direct 
consequence of this property, all the solutions of such systems are meromorphic functions of the complex 
time, and can be described explicitly in terms of theta-functions or generalized theta-functions. The property 
of meromorphicity led to the Kovalevskaya-Painleve integrability test, which was effectively applied to detect 
several integrable cases. In some cases, such as the famous Neumann system describing the motion of a point 
on a sphere with a quadratic potential, or the Steklov-Lyapunov integrable case of the Kirchhoff equations, 
the complex tori are Jacobians of hyperelliptic curves (or possibly coverings of the Jacobians). In other, 
more complicated situations, for example the Frahm-Manakov top on so(n) or the Kovalevskaya top, the 
vo. complexified tori are not Jacobians but Abelian subvarieties thereof, with a non-principal polarization (Prym 

\ varieties 

■ On the other hand, there are many other systems, including generalizations of the above ones, which 

preserve integrability but lose the meromorphicity of the complex solutions. In such systems the genus of 
the underlying algebraic curve (often the spectral curve) is greater than the dimension of the invariant tori, 
| and the latter are certain non- Abelian subvarieties (strata) of Jacobians. The algebro-geometric properties of 

such systems and the nature of the singularities of their complex solutions have been described in [Van995 , 
IAF001 IEPR031 IFG071 IEEKL9931 IEEKT994] . Until recently, however, the only examples known were of 
systems related to the strata of hyperelliptic Jacobians. 

In present paper we consider the first example of a mechanical system whose complex invariant varieties 
are strata of Jacobians of a non-hyperelliptic curve, here a trigonal curve of genus 3 given by the equation 
y 3 + p(x)y + q(x) = 0. The latter appear in the reduction to quadratures of the integrable Goryachev case 
of the Kirchhoff equations |Gor912] : the quadratures involve 2 points on the genus 3 curve (which has no 
additional involution in general, so not allowing reduction to a Prym) and so the quadratures lead to an 
incomplete Abel map which cannot be inverted in terms of meromorphic functions. This means that, as 
in other non-algebraic integrable systems, the Goryachev case cannot be detected by the Kovalevskaya- 
Painleve test. We emphasize that this example is not unique: it is in fact a member of a family of integrable 
Hamiltonian flows on the sphere S 2 that have extra cubic integrals and that were described recently in [VallOj, 
see also [VT12| . In particular this family also includes the non-trivial case found by Dullin and Matveev 
in [DM04] . As was shown in [Tsi051 IVT09 , most of the systems of the family are reduced to "trigonal" 
quadratures similar to those of the Goryachev systen0. For concreteness, the present paper considers only 
the Goryachev system. 



1 These can be related to other Jacobians in the case of dimension 2. 

2 There are several exceptions in the family: one of them is the classical Goryachev— Chaplygin system [Chap904 Gor916| . 
which is linearized on Jacobians of genus 2 hyperelliptic curve. 
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Our solution builds on the explicit description of Abelian functions of trigonal curves and (in the ter- 
minology of [BEL999 ) the more general (n, s)-curves. There has been a resurgence of interest in this area 
with many new analytic results obtained jBELOOl [BT051 IBG061 IEEMOP071 IEE091 INaklOl IEEQ1 lj including 
the inversion of the Abel-Jacobi map on strata of Jacobians [MP08, MP11, NYol2 . Although we will need 
to extend this work in various ways these new studies are foundational to providing explicit solutions for 
non-algebraically completely integrable systems. 

The paper is organized as follows. In Section 2 we reproduce the reduction to quadratures of the Goryachev 
system first made in |VT09i and interpret them as sums of two holomorphic differentials on a genus 3 trigonal 
curve C, also indicating its canonical form. The original variables of the system are then expressed in terms 
of coordinates of two points on the curve. In Section 3, following [EEMOP07 , we describe in detail the 
inversion of the complete Abel map on the (3,4)-trigonal curve. Here the main tools are the sigma-function 
of the curve and its logarithmic derivatives (Kleinian-Weierstrass functions) , which are direct generalizations 
of the Weierstrass elliptic a- and p- functions respectively (see [BEL9971 IBak897| ). We give a method of 
an effective calculation of the vector of Riemann constants for the trigonal Jacobians which allows us to 
calculate the sigma function explicitly by relating it to the corresponding theta-function. We also present an 
analytic description of the Wirtinger strata in the Jacobian of C as zeros of the sigma-function and of some 
of its derivatives. Section 4 describes the inversion of the incomplete Abel map. The formal explicit solution 
to the inversion problem is obtained from the formulae of the previous section (inversion of the complete 
map) by a certain limiting procedure. The resulting solution is given in terms of the sigma-function and its 
derivatives whose arguments are restricted to the 2-dimensional stratum in the Jacobian of C. Then global 
analytical properties of the solutions as functions of the complex time are described. We show that these 
functions have an infinite number of branch points and are single- valued only on an infinite ramified covering 
of the complex time plane. Finally the local singularities of the complex solutions are described in Section 5 
by using the expansion of the sigma-function near generic and special points of the curve C. An Appendix 
contains some rather long and technical proofs of two theorems. 



2. The Goryachev integrable case. Separation of variables and reduction to quadratures. 

Recall that the classical Kirchhoff equations describing the motion of a rigid body in an ideal fluid have 
(in an appropriate coordinate frame) the form 



where J = ( J\ , J 2 , and 7 = (71 , 72, 73) T are the angular and linear momentum respectively and H ( J, 7) 
is the Hamiltonian, which is also a first integral. In addition to the Hamiltonian the equations always possess 
the integrals (Casimir functions) 



Apart from the well known integrable cases of Kirchhoff, Clebsch, Steklov and Lyapunov (and their 
gyroscopic generalizations) there are some further special cases of integrability cases where an additional 
integral exists only under the condition C\ = 0. In the most classical case found by D. Goryachev |Gor912j 
and S. Chaplygin Chap904 the extra integral is cubic in J. In this case Chaplygin himself |Chap904| gave a 
separation of variables and reduced the system to quadratures containing integrals on a hyperelliptic genus 2 
curve. A detailed algebro-geometric description of the complex invariant manifolds was made in BvM987 . 
Further rather exotic special cases of integrability also exist for which neither separation of variables nor 
explicit solution were known until recently. Here we concentrate on the Goryachev case [Gor916 which was 
reduced to quadratures in |VT09] by using a bi-Hamiltonian structure and the corresponding separating 
Darboux-Nijenhuis variables. 




dH dH 



(1) 



Ci = (J, 7 >, <? 2 = <7,7>- 
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The Goryachev case, the focus of this paper, has Hamiltonian H\ = H and extra integral Hi that take 
the form 

Hi = J\ + jf + ~ Jf + — ^tt~' a ' & being arbitrary constants, 

3 7, 

(2) J3 

Hi = -\UJI + Jf + | Jf + + a4 /3 Ji. 

The corresponding Kirchhoff equations are then 

ii = | J2 h - ^72(071 + ^)73" 5/3 > 

J-2 = -| Ji J3 + a73 /3 + g7i(07i + & )73~ 5/3 > 

x 72 

J 3 = -a- 



(3) 7 3 2/3 ' 

4 

71 = 2(- J372 - J273), 

4 

72 = 2(--J 3 7i + J173), 

73 = 2(71 Ji - 72 Ji). 

Without loss of generality one can set Ci — (7,7) = 1. Then, since C\ — 0, equations ([I]) can be reduced 
to a Hamiltonian system on the cotangent bundle of the unit sphere S 2 = {(7, 7) = 1} with coordinates and 
momenta 

J 1^1 — J271 

(4) " = 73, Pu = 2 2 — , 4> = arctan(7i/7 2 ), p<p = -.h- 

7i + 72 



In terms of these the original variables become 



n sin < > iu + cos 0(1 — u 2 )p u u cos cj)p^ — sin 0(1 — u 2 )p v 



Jl — n > ^2 — 



u 



2 



(5) xi = yl — u 2 sin0, xi = yl — u 2 cos </>, 23 = u . 

The paper |VT09j introduced separating variables q± , qi as the roots of the polynomial 



A(A) = A 2 +^_^_^ A _ 



uie 



Observe this polynomial depends not only on the coordinates u, <f> of S 2 but also on their momenta p n 
Following |VT09] . in the Darboux coordinates qi,Pi on T* S 2 such that 

(6) UhPj} = 4/> {ft, Qj} = {Pi,P]} = 0, 

one has 

3/2 



P2 - Pi 
P4> = «9l92 : 

(7) qi - q2 



2 gipi - q 2P i 

3 91-92 



o '/'» ) = -9.1. - 92- ' 7= = „ = -9192, 



,1/3 f 



and the two Hamiltonians take the Stackel form 

Hi = (5- 1 ) 11 [/(p 1 ,9i) + (S-^nUip^qi), 
(8) #2 = (S-^uUipi, qi) + {S- l )iiU{p 2 , 92), 

with the Stackel matrix 

1 1 



S = 



3i/(2pi9i) 3i/(2p 2 q 2 ) 
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Setting for convenience A^ = qi, fJ-i — 2/3iqiPi, we get 

U(q,p)^U(\^) = ^(X^ + X 4 + a 2 /4-bX 2 ), S=(^^ 

and then observe that the above relations are equivalent to 

1 



H 1 -H 2 — = U(\ llt i 1 ), Hi~H 2 — = U{X 2 ^ 2 ) 
Mi M2 



Then, upon fixing the values of the integrals, Hi = hi, H 2 = h 2 , the pairs (Xi, Mi) are subject to the algebraic 
relation 

/j,hi - h 2 



+ X 3 + -— - bX 
4A 



or 
(9) 



F = X 4 -bX 2 + (m 3 - hifji + h 2 )X + a 2 /4 = 0. 



As was also mentioned in |VT09j . equation ([9]) defines an algebraic curve C C C 2 (A, //), which, for generic 
values of hi,h 2 is smooth, has genus 3, and is non-hyperelliptic, i.e., cannot be transformed to the form 



y 



P(x) by a birational change of coordinates. A basis of holomorphic differentials on C is 



(10) 



dX XdX fidX OF 2 



2.1. The quadratures. Let ti,t 2 denote the time of the flows on T*S 2 defined respectively by the Hamil- 
tonians Hi and H 2 . To describe the evolution of qi = Xi with H 2 , we use the bracket ([6]) and the expressions 
J5J to obtain 



— Xi 

dt 2 dpi 3 



dH 2 8H 2 2 , 

— 1X1 = 



/'2 



MlM2 



\Ml~M2 (M1-M2) 5 



— jAi 

3 Mi - M2 



(C/i - C7 2 ) + 



Ui + Mi^— (mi^i - M2E/2) 

d/ii ^1 - [i 2 



M1M2 dUi 
Mi - M2 9mi 

t/^^A^Mi) 



2 . 
3 ZAl 



ffi(A.M) 



/ 



In view of the expression for F in ^ this is equivalent to 
(11) 



dXi 2 M2 dF(Ai,Mi) , . , dX 2 2 mi <9F(A 2 ,m 2 ) 

— — = —1 , and, similarly, — — = —1 . 

dt 2 3 Mi _ M2 »Mi dt 2 3 M2 - Mi 0M2 



Similarly one obtains for the flow with the quadratic Hamiltonian Hi 
dXi 2 1 <9F(Ai,mi) dX 2 2 1 



(12) 
and also 
(13) 



9F(A 2 ,m 2 ) 



dti 



3 Mi - M2 c^Mi 



dti 



3 M2 - Mi ^M2 



rfMi _ dfiidXi _ 2^ 1 <9F(Ai,mi) d[i 2 _ 2^ 1 dF(X 2 ,^ 2 ) 



dti dXi dti 3 mi — M2 <9Ai c?ii 3 M2 — Mi 

Expressions (fl"2j) . (fTTj) yield the following quadratures in the differential form 



9A 2 



(14) 



dXi dX 2 2 

W{X^{)Jd^i + 9F(A 2 ,M 2 )/9M2 ~~3* * 2 



Mi<^A 2 



H 2 dX 2 



dF(Xi^i)/dfii dF(X 2 ^ 2 )/dfi 2 



-1 dti. 



We will return to these later. 
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2.2. The original variables in terms of the separating ones. From ([7]) and other formulae from the 
paper |VT09j one has 

2 _ _( Mi - M2 \ 3 2/3 _ Mi - Mg 

^3 - t r~ ' ^3 



^3 = 



Ai — A 2 / ' Ai — A2 

3 Ai^ 2 - A2M1 



2 Ai - A 



2 



2t. , / 1 , / A*i - M2 N 3 



exp^) = — AiA 2 Wl , , 

a y \Ai — A; 

72 +«7i = — A1A2 ( 1+ ( ^ 7^) I, 72-«7i 



a \ \Ai — A 2 / / ' " 2«AiA2 



(15) 



Ji + «J 2 = -« f - - — ^2?V>J e v^l~ 
a Ai + A 2 / Mi _ M2 



2 A1A2 \ Ai — A2 

Ji — zJ2 = ((Ji + «J 2 )(72 + «7i) + 2iJ373) (due to the condition J • 7 = 0) 

72 - »7i 

2AiA 2 (fl + ^^C) (A 1+ A2)+ 3( ^ 1 - At2)2(AlAl2 -^ A2)N 



(Ai - A 2 )V v z; (Ai - A 2 ) 3 



Ml— M2 

Ai -K2 

Substituting the above expressions into the integrals @ and using the equation of the curve (O for each 
pair (Xi,fii) one may verify the identities H\ = hi, H% — h 2 . Thus to solve for the original system ([I]), ([2]) it 
suffices to solve for the pairs (Ai,/ij). 

2.3. Canonical form of the curve and of the quadratures. Rather than using the variables (Ai,/i;) 
directly we now make one final birational transformation to bring the curve C to a canonical form. This 
allows us to make connection with the literature on (n, s) curves and so permits us to solve for the motion 
in terms of the multi-dimensional cr-function. 
By making the birational change 

the curve ([9]) can be transformed to the canonical trigonal form with respect to y 



hi 2 (4,, 2y/2h 2 3 b 2 
— x y—\x-\ 3-^— x - 2-x + 1 



(17) G{x,y) = y A -2^-x I y~ ( x 4 + ~^x A - 2-x 1 + 1 ) = 0. 



a 



We will refer to this curve as C. In the terminology of [BEL00 (see also the next section) this is a (3,4)-curve 
having one infinite branch point 00, where all the 3 sheets of the covering C — > P 1 = {x} come together. 
Under the above transformation the holomorphic differentials uj\ , 012, W3 in (|10|) take the following respective 
forms, 

xdx dx ydx dG 2 2h x 2 

(18) n ^-dGWv ni = dGld~y n * = dGld~y ~dy- =3y -— x - 

Then, in the new coordinates, the quadratures (|14l) take the form 

x\ dx\ x 2 dx 2 2 



idt 



2. 



dG(x 1 ,y 1 )/dyi dG{x 2 ,y 2 )/dy 2 3 
dG{x 1 ,y 1 )/dy 1 dG{x 2 ,y 2 )/dy 2 3 ' 
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2.4. Non-algebraic integrability. We are now in the situation described in the introduction. The genus 
of the curve C is greater than that of the dimension of generic invariant tori of the system: here we have 
2 separating variables while C is of genus 3. Such a situation occurs in many algebraically integrable 
systems, for example the Clebsch integrable case of the Kirchhoff equations or the Kovalevskaya top. In 
these examples however, although the genus of the underlying curve (the spectral curve of the corresponding 
Lax representation) is greater than the number of degrees of freedom, the relevant curves possess an additional 
involution which extends to the Jacobian variety. Then, as a rule, the complex invariant manifolds of the 
systems turn out to be 2-dimensional Abelian (Prym) subvarieties of the Jacobians, whose real part gives 
the invariant tori. This however is not the case for the Goryachev system on T*S 2 we are considering. For 
generic hi, h 2 the curve C has no further symmetries, and the differentials in (fl~9| do not reduce to those of 
a genus 2 curve. As we shall see below, this pathological property means such systems are not algebraically 
completely integrable. In particular, their complex invariant manifolds are non-Abelian subvarieties of the 
Jacobian of C and the variables Ji, 7$ are not meromorphic functions of the complex times t\, t 2 . 

In terms of the differentials (fT5|) the quadratures (|14l) may be extended and written in the following 
integral form 

u 2 = -2it 2 /3 + u 20 , p t = ( XhVi ) g c , 
u 3 = -2iti/3 + u 30 

where 1620,^30 are constant phases and the coordinate u\ is a transcendental function of 1*2,^3, whose 
properties will be described in the next sections. Thus (|20[) defines a map of the symmetric product C x C 
to a codimension one subvariety (stratum) of the Jacobian variety of C. To invert the map, i.e., to express 
symmetric functions of the coordinates (xi, yi), (x 2 , y 2 ) and consequently the variables J^, 7, in terms of 
u 2 , U3, at least locally, we next recall some basic facts about the standard Jacobi inversion problem associated 
to trigonal curves. 



3. Jacobi inversion problem for the genus three trigonal curve 

The curve (fT?]) belongs to a class of (n, s)-curves. These are smooth curves with s > n > 2 and gcd(n, s) — 
1 that have one point at infinity and whose affinc part may be defined defined by an equation 

V n ~ x ° ~ v an +p s x a y P = 0, 0<a<s-l, O<0<n-1. 

a,/3 

The curves are of genus g = (n — l)(s — l)/2 and their properties and relation to integrable hierarchies of 
KP type are widely discussed in the literature. 

Below we concentrate on the trigonal (3,4)-curve C, which we write in the canonical form 

(21) f(x, y) =y 3 + ([i 2 X 2 + fl 5 X + [l 8 )y - (x 4 + H3X 3 + H 6 X 2 + H 9 X + ^12) = 0, 

where [ij are parameters. This is more general than our curves and we will specialize in due course. 
The curve has one infinite point 00, where all 3 sheets of the covering C — > P = {x} come together. Let 
£ = a^ 1 / 3 be a local coordinate in a neighborhood of 00. That is, £(00) = 0, and the coordinates in the 
vicinity of this point have expansions 

Choose the vector of holomorphic differentials ft = (f2i, £l 2 , ^3) T , 

(23) n = 
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Near oo, they admit the expansions 

^i = -e 4 (i + fe 2 -^ 3 + o(c 5 )))^, 

(24) n 2 = -e (i + f e 2 -^e + o(e 5 )) 

^ = -(i-fe 3 -fe 4 + o(e 6 ))de, 

so the orders of their zeros at oo decrease. 

Next, choose a canonical basis of cycles of H\{C, Z) 

(o, b) = (ai, a 2 , o 3 ; bi, b 2 , b 3 ), such that o; o ctj = bi o b, = 0, Oi o bj = -b» o aj = 1 

and introduce matrices of periods of the above differentials 

(25) A=(<ffk) , B=[l SlA 

\ Ja J /t,i=i,2,3 \ Jb > l .... :.J .i 

Throughout the whole paper we will use two normalization of periods: the first, given above, [A, B), 
where we have specified the differentials; the second utilizes the so called a-normalized differentials for which 
the matrix of periods takes the form (I3, r). Here r is the Riemann period matrix 

(26) t = A~ x B, t t =t, lmr>0. 
We denote the a-normalized holomorphic differentials by 

(27) Tl = A~ 1 n & f ti t = 5ij, i,j = 1,2,3. 

The first normalization is used in the definition of the c-function while the second in the definition of 
Riemann's theta-function 0(z;t) associated with C and the Riemann period matrix r given in (|26|) : it is 
defined by the series 

(28) 9(v; t) = exp {tirn T ■ r • n + 2nrv T • n} . 

With these normalizations we define the Jacobi variety Jac(C) of the curve C to be the quotient Jac(C) = 
C 3 (wi, u 2 , u 3 )/{A © B} and also Jac(C) = C 3 («i, v 2l v 3 )/{l 3 r}, where the vectors u and v are related by 
v = (v 1 ,v 2 ,v 3 ) T = A^u. 

For a positive divisor T> = P\ + P 2 + P3, Pi = (xi,yi) £ C, we consider the Abel map with a base point 
P ,W:C xC xC — > Jac(C), given by 

/•Pi [P2 fP 3 [Pi _ [P* _ rP S _ 

(29) 21 : X> — > n + / fl + / n = u, or, equivalently, / ft + ft + fl = v, 

J Pa J P a J Pa J Pa J Po J Pa 

Henceforth we assume that Po = (00, °°)- The Jacobi inversion problem refers to the inversion of Abel map. 
For the 3 = 3 case being considered this is solved in terms of Riemann's 0-functions as follows. 

Theorem 3.1. (Riemann) Let e be a vector such that the junction F{P) = 8(f^£l — e;r) does not vanish 
identically. Then F(P) has exactly 3 zeros on C, P\ 1 P 2l P3, and these provide a solution of the Jacobi 
inversion problem 

fPl+P2+P 3 _ 

(30) / fl = e Koo, 

J 3oo 

where Koo = (Ki, Ki, K3) is the vector of Riemann constants with base point 00 and whose components are 
given by 



-)-^g(£n, ( p)/>), 



(311 K, = i(2^ + r rr ) - — ^ ( / O/i/'i / <),. I . ,- = 1.2.3. 

l^r 
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The corresponding divisor T> = Pi + P2 + P3 is non-special, and in the vicinity of T> the map 21 is uniquely 
invertible. 

The rank of the Abel map is maximal on non-special divisors and decreases on special subvarieties of the 
Jacobian, the Wirtinger strata. Here W^°' = C C C Jac(C) where 



(32) 



rP_ 

W {1) : veJac(X), v= n + Kcc, VP 6 C 

J OO 



oc 

W {2) : veJac(X), v= / « + / O + K^, V(P X , P 2 ) € C X C . 

00 j 00 



The equation 9(v;t) = defines a codimension one subvariety € Jac(C) (with singularities for g > 2) 
called the theta- divisor, which coincides with the stratum W^ 2 \ This is equivalent to the fact that for any 
points Pi, P 2 G C, 




(33) / « + / O + Koojr =0 



Note that a consequence of Riemann's theorem is that the vector K oo for which (|33)l holds is unique. 
We also introduce characteristics, represented by real 2x3 matrices 

£ T \ _ ( £1 £2 £3 
e' T )-{e[ £ 2 £ 3 

so that any vector u g Jac(C) can be written in the form u = Ae' + Be, e, e' € R 3 . We denote the 
characteristic of u by [u]. In what follows we concentrate on rational characteristics, in particular half- 
integer ones, for which £j,£j = \ or 0. We shall also need the Riemann theta- functions with characteristics 
\e\ given by 



£ T 



(«; r) = exp {«tt(£ T t£ + 2£ T (i; + e')))} 6>(u + r£ + e'; r) 

= exp {i?r(n + e) T r(n + e) + 2m(n + e) T {v + e')} 



3.1. Calculation of Koo. The vector of Riemann constant K oo given by ([31]) includes Abelian integrals 
over o-cycles, which are difficult to calculate directly. It is known that if a curve has a point P* such that 
the canonical divisor is linearly equivalent to 2(g — 1)P* then the vector of Riemann constants with base 
point P* is a half-period |Fay973| |FK980| . Such is the case for a hyperelliptic curve. For all {n, s)-curves 
Nakayashiki [NaklO observed that these admit a holomorphic differential that vanishes to order 2g — 2 at 
the point P = (00, 00) and so we have the following. 

Proposition 3.2. The vector of the trigonal curve C is a half-period in Jac(C). 

We may further restrict the choice of K^. There are 64 half-periods in Jac(C), 28 odd and 36 even ones. 
Up to exponential terms the a- function on C is the 0-function shifted by -ftToo (see below) . Thus the leading 
terms of the expansion of the 0-function at K^, coincides with those of the corresponding er-function. For 
(3, 4) curves it was shown in [BEL999] that the cr-function begins with an odd order Schur function and it 
follows then that Kqo is an odd half-period. Of course the explicit expression for depends on the choice 
of homology basis on C and so to proceed further we must first fix the homology basis. We may specify which 
half-period then corresponds to Koo by first calculating the period matrix and then checking the condition 
(f33|) by computing the expansion of 0-function near all odd half-integer characteristics. For example, for the 
curve C given by 

y 3 - Ayx 2 - a; 4 - 3a; 3 - Ax 2 - x + 5 = 
and in the Tretkoff-Trctkoff basis of cycles on it given by algcurves package of Maple one obtains 

1/1 1 



(34) [*~]= 2 l 1 1 
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We note in passing that we already encounter here one of the differences with Jacobians of hyperelliptic 
curves: for a genus three hyperelliptic curve the vector of Riemann constants is given by an even singular 
half-integer characteristic whereas the characteristic (|34[) is non-singular and odd. 

3.2. The er-functions. Apart from the theta-functions, in many cases it is more convenient to use the a- 
function of u. To describe it, we first introduce the basis of 3 meromorphic differentials Y = (Ti, T2, T3) 7, 
having a unique pole at the infinity of C and specified by the pairing conditions 

(35) Resp =00 (U(P) [ Tj= Sy, i,j = 1,2,3. 
Then, following [EEMOP07] . 

dx I Tl 

(36) T = — ^ -2xy + ^x 2 
where 

Ti = -(5x 2 + (/Z1/Z2 - 3^3)2; + ^2^4 + Ve)y + mu 2 y 2 + 3^ix 3 

2 , n r, \ 2 / 1 \ . ^ , . / 1 1 n 



- (ji 2 + 2(i 3 m - 2fi 4: )x - (fi 5fJl2 + fiQ^ + fi 3f i 4 )x + -fiif x (x, y) - ^ -(is - -Mi ) f y {x, y). 
The matrices of periods of the differentials 



(37) S = - If dTi , T=-U dT, 

\Ja k / i,k=l,2,3 \Jb k / i,fc=l,2,3 

satisfy the generalized Legendre relation 

Let us also introduce the normalized second period matrix x = SA^ 1 , which is necessarily symmetric. 
The fundamental cr-function of the curve C is defined by the formula (see e.g., [BEL997, EELOO ) 



(39) 



where 6[K](v) is the theta- function with the characteristic corresponding to the vector of Riemann constants 
(for a chosen base point and homology basis), and C is a constant depending on the period matrix A and 
the coefficients of the curve. The constant C provides the modular invariance of (|3T)1) . just as in the case 
of the Weierstrass elliptic a-function. An explicit expression for C is given in [EEMOP07 and it is not 
necessary for our exposition (we will deal only with ratios of the sigma- function derivatives). According 
to the definition, the fundamental cr-function is normalized in such the way that its expansion at u ~ 
starts with with the Schur- Weierstrass polynomial (see BEL999 for details). It follows that a(u) is just 
the theta-function of C whose rescaled argument is shifted by the vector Koa and multiplied by a quadratic 
exponent of u. Thus, the knowledge of the corresponding characteristic [K\ calculated above is important 
in the explicit description of a(u). 

The function o(u) is quasi-periodic just as is 6(v): when u is shifted by a period An + Bm, n,m e Z 3 
then a(u) is multiplied by the exponent 



(40) o(u + An + Bm) = cr(tt)cxp j (Sn + T m) T \ u + -An + ^Bm 

For the case of the curve (f2Tj) , at the origin u = 0, the cr-function admits the following expansion (see 
[EEMOP07p 

(41) cr(tt) = u x - u%u 3 + ^rul - J^ u l + y u 3 u 2 + ^§ U 3 U 2 - y"3 w 2 + higher order terms. 



1 . 1- 
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It follows from property ([55)1 and the definition (HE)]) that for any points P\,P 2 G C 
(42) erf / fi+ / f2]=0. 



That is, in M-coordinates on Jac(C), the stratum W^ 2 ' is given by the condition a(u) = 0. 

3.3. Inversion of the Abel map in terms of the sigma function. We next introduce the Kleinian 
multi-index symbols 



These are multiply periodic (Abelian) functions 

(44) pj(u + An + Bm)=pj(u), 

where J is arbitrary multi-index with more than one entry. It will also be convenient to denote throughout 

OUi 

though these are not Abelian functions. 

Then, as was first shown in BELOO, EEL00 (though in a somewhat hidden fashion), the problem of 
inversion of the map (|29[) is reduced to solving the following two equations with respect to x and y, 

(45) p 33 {u)y + p 23 (u)x + p 13 (u) = x 2 , 

(46) (P23(«) - P333{u))y + (p22(«) ~ p233(u))x + p 12 {u) - pl 33 {u) = 2xy, 

the solutions of which give the coordinates of the points Pi,P 2 ,P 3 in (f2T?l) . In contrast to the inversion 
problem for hyperelliptic case, equations (|45[) and (|46p both contain the variables x and y. By elimination 
of y, one obtains a cubic equation with respect to x, whereas elimination of x yields a cubic equation for y. 
In the first case we get (the new expression, though a simple consequence of the previous) 

V(X]U) = 2x 3 - (3p23(«) - p333{u))x 2 
^ - (P33(u)p22(u) + p23(u)p333(u) - p33(")p233(«) ~ P23(") 2 + 2p 13 (u))x 

- P33(u)(pl2{u) - p 133 (u)) + p 13 (u)(p 23 (u) - p 333 (u)) 

= 2(x - xi)(x - x 2 )(x - x 3 ). 

These formulae enable one to express the elementary symmetric functions of Xi as Abelian functions of u. 
We shall use some of them in the next section. 

3.4. The inverse trace formula. Let c/, = (0, j/j), i = 1, 2, 3 be the three points on the curve C over x = 
and Ti = f q ' ft be their images in Jac(C). By applying the classical method of residues one can also obtain 
the following "inverse trace formula" : 

111 { 1 ' ' 

(48) — + — + — = Y log a(r t - u) + k, V, = I — r, 0, 



xi x 2 x 3 ^ ' \f y (0,y l )' ' f v (Q,yi) 



Here Vi is the tangent vector to C Jac(C) at the point r<, and k is a constant depending on the curve 
C only. 

Indeed, consider the single-valued function F(P) = 9 ( ft — v — Ko^j defined on a simply connected 
dissection C of C with boundary dC. Then, for a mcromorphic function f(P), P £ C with the poles Qi, ■ ■ ■ ,Q S 
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on C, the residue formula gives (see, e.g., |Dub81 . BBEIM94]) 

(49) /(Pi)+/(P 2 )+/(P 3 ) = ^- <f f(P)dlagF(P) -J2Re SQ J(P)dlogF(P). 

Set here f(P) = 1/x and observe that this function has simple poles precisely at Qj = i = 1, 2, 3 and that 
the expression is independent of the choice of local parameter. Then, using the expansions of F{P) in the 
neighborhood of and the relation (f39|) between 9(v) and cr(it), one arrives at (|48l) . Alternatively, one can 
derive a cr- function analogue of formula (|49|) . 

Notice that Vi + V2 + V3 = 0, which, in view of the quasi-periodic property of a(u), ensures that the right 
hand side of (|4"5|) remains unchanged when u changes by a period vector of Jac(C). The constant k can be 
calculated explicitly by letting x\,x 2 ,X3 — » 00 in (|4"51) and evaluating the right hand side (see also formula 
d5l| below). 

For the curve relevant to the Goryachev case under consideration we have simplifications. 
Proposition 3.3. In the special case p% = 0, p\ 2 — 1 in the equation of the trigonal curve (|2ip one has 

(50) Q t = (0,p l ), P = exp(27rz/3), V; = i 0, 1) = i(p l , 0, p 2 *). 

In i/ws case formula (|48[) tafces i/ie /orm 
1 1 1 1 3 

(51) — + — + — = - V (p^i log tr(r< - «) + p 2i <9 3 log ^(r, - u)) + fc, 

Xl X 2 £3 3 f-f 

z— 1 

(52) . = ^(^ + p 2 ^ + p 4 ^V ri =r pi n. 

3 V 0-2(7-3) cr 2 (ri) 0-2 (r 2 ) / 

Proof of Proposition 13.31 The formula (|5~Tj) follows from (|48l) under the conditions (|50"|) . To calculate the 
constant fc, let the points Pi, P 2 , P3 € C to tend to 00 in such a way that under the Abel map (|3"0"|) one has 
ui = U3 = 0, and the coordinate u 2 tends to zero. Then the left hand of (1511) tends to zero, while the right 
hand side becomes the sum of k and of the terms 

p 2l oi{ri) + pV 3 (n) 
3o-(rj) 

Now as <7i(r,) = 03(1^) = o~(ri) = this expression in indeterminate, but applying l'Hopital's rule to each 
term we obtain 

p 2l (Xi(ri + u 2 e 2 ) + p i <ra(n + u 2 e 2 ) .. p 2l o- 12 (ri + u 2 e 2 ) + pV^r; + u 2 e 2 ) 
una ; ; = lim 



u 2 ->o 3CT(r 3 + u 2 e2) "2^0 3cr 2 (ri + u 2 e 2 ) 

p 2 Vi 2 (r 4 ) + p i a 23 (r i ) 



3a 2 (ri) 

where e 2 = (0, 1, 0) T . The latter fractions are well determined. Finally, as we will prove in Proposition 4.2 
below, o-\ 2 (ri) = for any i = 1, 2, 3, and so ([53)1 gives (|52")) . □ 

4. Solving the inversion problem on the stratum and sigma-function solutions of the 

Goryachev system 

We now identify the trigonal curve (|TT)) that appears in the quadratures for the Goryachev system and 
the curve C in (f2"Tj) by setting 
(54) 



pi = /i 4 = 0, p 2 = -2hi/a, /i 5 = /i 8 =0, p 3 = 4/i 2 I-) , p 6 = -2b/a, p 9 = 0, p,\ 2 = \. 

Notice again that the map C x C — > Jac(C) in ([2"u]l contains only 2 points on C and cannot be identified with 
the full Abel map (|2T))) . According to the previous observations, (|20p maps the symmetric product C x C to 
the stratum C Jac(C) given analytically by the condition a(u) = 0. 
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4.1. The inversion of the Abel map on the stratum One can extend this map by adding a third 

fixed point on C, in particular, one of the three points qi = (0, p l ) G C: 



P1+P2 



2oc 



f2 = U + Ti 



a, 



u being the right hand side of (|20l) . Then, by using (|4"T1) , one obtains the following formal complex expressions 
for the symmetric functions of x\, x 2 (which hold for any i = 1, 2, 3 on the right hand side): 

(55) Xi + X 2 = 3p 23 (u + - p 333 {u + Ti), 

X\X 2 = p 33 (u + ri)p22(u + Ti) + p 23 (u + ri)p 333 (u + n) 

(56) - p 33 (u + r l )p 233 (M + r-i) - pl 3 (u + n) + 2p 13 (u + n), 

where, according to (j2"0")l . u 3 — —2itx/3, u 2 = — 2tt 2 /3, and the coordinate u\ is defined (but not uniquely!) 
from the transcendental condition cr(u) = 0. 

On the other hand, letting x 3 — > oo in the inverse trace formula (|5ip . we hnd 

1 3 



X\ + x 2 



X\X 2 



r(«)=0 



1 3 



P 2 '0-2,3(T-») 
0-2 (U) 



Upon combining the above with (|55|) . one obtains the following alternative to ([5B]) 



(57) 



3p23(« + - p333(« + n) 



a(u)=0 



V. t = -(p\0,p 2 *). 



Ei=i ^ log o-(m - n) + fc 

There is another way of writing the formal solution to the inversion of (|20|) . (For the case of hyperelliptic 
curves this was proposed in |Gra990| IJor992j and also used in |EPR03| .) Namely, consider again the full 
Abel map (|29|) with the three points Pi — (xi,yi),P 2 — (x 2 ,y 2 ),P 3 = {x 3 ,y 3 ) and observe that 

XiX 2 + X 2 X 3 + X 3 Xi 



X\ + x 2 = lim 

X3— >oo 



Then, in view of ([47]). 



X\ + x 2 = lim 

P 3 ->oo 



X\ + £2 + 2:3 

P330")P22(m) + P23(")p333(«) ~ P33(«)p233(w) ~ P23(") 2 + 2pi 3 (u) 



3p23(«) - p333(«) 

Using the definition of pij in (1431) and taking the limit (for which a(u) — 0) one obtains 



1 <T 3 (m)<7 2 , 3 (") " CT2(w)g3,3 (m) , 1 g2(M)03(M) ~ &3(u)(t) 2 (u) 

2^ U) al(u) + 2 a 3 » 



(58) X\ + x 2 
where we set 

(59) 4>{u) = <j 2 (u) - ct 3 , 3 (m), ^•(u) = £7-j )2 («)-o-j,3,3(«) ) j = 1, 2, 3, ueW {2) , 



and m is given by the right hand side of (1201) . In a similar way we find 
(60) X\X 2 



1 ^ (Ti(m)o-3.3(m) - Q- 3 (M)gi. 3 (M) | 1 O-3(m)0i(m) - CTi(m)03(m) 

2 0W oS(«) + 2 of(«) 



it G 



IU< 2 



Next, using the relation (|45[) for the pairs (a?x,j/i), (0:2, 2/2)1 we find 
(61) y-j = (x 2 j - p 23 (u)x, - p 13 (u))/p 33 (u), 

and, therefore, 

x 2 yi - xiy 2 a:ia;2 - pizju) yi - y 2 _ Xi + x 2 - p 23 {u) 

x 2 -xi P33(") ' X 2 -Xi P33(") 

where u is again given by the right hand side of ([20|) . In view of (|43|) and the condition a(u) = 0, this gives 
the following new compact expressions 



(62) 



x 2 yi - xiy 2 
x 2 - x\ 



Q"i(m) 
03 (m) 



cr(li)=0 



2/1 - 2/2 

X 2 - Xl 



Q-2 («) 

03 (m) 



cr(li)=0 
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The above formulae lead to the following analytic description of the Wirtinger strata defined in (j3"2")l . 
Proposition 4.1. (1) The strata W {Q \W {1 \W {2 ^ are given by the conditions 

W {0) : a{u) = a 3 {u) = a 2 {u) = 0, 
(63) W {1) : {u\ a(u) = 0, a 3 (u) = 0} , 

W {2) : {u\ a(u) = 0} 
(2) The coordinates x,y of the curve C — C Jac(C) admit the parameterization 



(64) xi =-^T, yi = g i( tt K3(t»)-^(t») gl , 3 (u) u€W p } 



Proof. The description of was already given by (|42l) . To pass to the stratum W^ 1 ' one should let 

x 2 — > oo. Then both sides of (l58l) and (|60| tend to infinity, which happens if and only if a(u) = and 
a 3 (u)=0. 

Next, in view of (|58j). (|60|). and the condition 0-3(14) = 0, 

,. ^1^2 0-1 (u) /" Pl 
hm = 2:1 = 7 — r, u= it, 

x 2 ^ooX 1 +X 2 0- 2 {U) 

which gives the first expression in (|64)) . Letting here 2:1 — > 00 gives a 2 (u) = for u = 0. The latter also 

follows directly from the expansion ([41"T) . The second formula (|64| is obtained by the appropriate limit from 
(EH). □ 

We shall also use the following 

Proposition 4.2. Let w = u — r$ = (w\, w 2 , w 3 ) and set k = o~j k(rj). Then 

cr| — a 3 — and the first few coefficients of the expansion near Vi € i = 1,2,3 

/ \ (i) 1 (i) 9 («) 1 U) 1 1 (i) 2 («) fil 

cr(rj + w) = o\ w 2 + — cr\ \ W\ + a\ 3 W\ w 3 + — cj 3 3 + — a\ 2 + °i 2 w i w 2 + °2 3 w z w 3 

(65) +- tr|*\ ) ! wi 3 + - er^ 2 iu 2 w x 2 + - er^ 3 w 3 -u^ 2 + - er^ 2 i0i w 2 + a[^ 2< 3 w 3 Wi w 2 

1 (i) 9 1 (i) 3 1 fi") 9 1 fil 9 1 (i) 3 j, 

+ 2 °1, 3,3^3^1 + £02, 2,2^2 + 3 CT 2, 2, 3 W 3 ^ 2 + - CT^ ^ 3 W 3 W 2 + - (T^ 3 + 0(10 ) 

are related as follows 

(66) 4 3 = 4\ ^3 = -pV«, ^3,3=341 

(67) ^=^1,1=0, <Ti:i, 3 = ^i 3 = -2pvg. 

Here, as above, p — exp(2-7rz/3). 

The proof uses the expansions of the coordinate y near the points qi = (0,p J ) G C, 

2/ = P l + \^2{p l + 1) + ^P l )x 2 + ^pV 3 + 0(a; 4 ) 

(68) 3 3 

= p< + -{-p 2P 21 + M6^)^ 2 + tfp'x 3 + 0(x 4 ). 

and the corresponding expansions of the holomorphic differentials (p~8|) . Thus, for example, near (0,p) one 
has 

n l = U- l - 2p2 ^-^ x 2 -\ P p^ + ^ p2 ^- 2 ^^ x^ o&)) dx, 

\6 9 p 9 27 p / 

(69) 2 = (f x - I , 3 - ^ , 3 + 1 (3pV6 p V2)M6 J + O^)) dx, 

' P 2 1 2 2 1 2 3 , 1 P 4 ^6 2 +P 2 P6H2 - P2 2 4 , ^ V 
P UfiX /i3/3 x H X + C X 

3 9 P p 9 p ^ 27 p 2 V 

and the expansions near (0,p 2 ), (0, 1) are obtained by replacing p above by p 2 and 1 respectively. 
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The remainder of the proof is technical and is presented in the Appendix. 



4.2. Analytic properties of the inversion on W^ 2 K One should stress that formulae ([58|) - (|62l provide 
only a local analytic solution to the inversion of the incomplete Abel map ([20]) . This is because cr(u) = 
is a transcendental equation and the argument ui is an infinitely- valued complex function of u 2 ,U3. This 
fact also admits a geometric description. Namely, let W^ 2 * 1 C C 3 (iti, 1*2,1*3) be the universal covering of the 
stratum C Jac(C) and consider the projection 7r : W^ 2 ) — > C 2 (u2,u^). 



Proposition 4.3. Assume that Jac(C) has no Abelian subvarieties. Then, under the projection it, the variety 
W^ 2 ) is an infinitely-sheeted covering 0/ C 2 (ii2, U3) ramified along the subvariety W C W^ 2 ' defined by the 
conditions {cr{u) = 0, 01 (it) = 0}. 

Moreover, let {u*} C W^ 2 ' be the equivalence class corresponding to any point u * G W^K Then the 
projection tt{u*} forms a dense set on <C 2 (u 2 ,U3). 

Note that in the case of hyperelliptic curves of genus 2, when the codimension one stratum coin- 
cides with the curve itself, a similar description was made by Jacobi in connection with the inversion of a 
single hyperelliptic integral (see Mark992 ) , whereas for hyperelliptic curves of any genus and strata of any 
codimension a similar theorem was proven in |AF00| . 



Proof of Proposition \4-3\ The proof follows the same lines as that in [A~F00 . Namely, let t>i,...,t)g G 
C 2 (ti2,ti3) be the 7r-projections of six independent period vectors of Jac(C). For any point u* G W^ 2 \ the 
projections of its equivalence class on W^ 2 -* has the form 



|tt(m*) +YjmjX>j I nij G zj . 



Under the condition of the proposition the periods themselves are not commensurable and consequently the 
integer coefficients mj can always be chosen in such a way that the above sum will give a point in any small 
neighborhood of any point of C 2 (1(2,1*3) fixed a priori. Next, since W^ 2 -* is defined by the transcendental 
equation o~(u) = 0, for any point (1X2,1*3) G C 2 there exists an infinite number of solutions u\. In other 
words, the covering 7r : W^ 2 -* — > C 2 (u 2 , U3) has an infinite number of sheets. Finally, by the implicit function 
theorem, 7r is ramified over the points satisfying (71(11) = 0. This establishes the proposition. □ 



4.3. Sigma-function solutions of the Goryachev system. We conclude the section with the formal 
cr-function solutions for the original variables of the Goryachev system. First, applying the transformation 
(|TB| to ([13, one gets 

2 (x 2 yi-xiy 2 \ i 3 [a y 1 - y 2 
7 3 = > j3= o\/o ' 

\ X 2 ~ X\ J I V 2 X 2 — X\ 

1 / ( x 2 y\ — Xiy 2 \ 3 \ 

72 + «7i = 1 - , 72 - 171 = -1x1X2, 

x\x 2 \ \ x 2 - Xi J J 

1/2 

a, , , x 2 yi - x x y 2 



(70) 



Jl + iJ 2 =-\ ~{xi + x 2 ) 

V 2 \ x 2 - xi 

Ji - iJ 2 = (( Ji + iJ 2 )(72 + «7i) + 2iJ 3 73) ■ 

72 - 171 
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Then, upon comparing the above with the sigma-function expressions ([5 
simplification 



(ISOT) . ([52), one obtains after 



(71) 



where we set 



2/3 
7 3 



72 - ^ll 

72 + 111 

Jl + tj 2 



*i(«) 
$ 2 (w) 



CTi(u) 



0"3 O) 
I $ 2 (lt) 



J 3 = ^ 7 



cr(u)=0 



a (7^1/. ) 



2 V 2cr 3 (tt) 



r(u)=0 



2i 



2 ^ 3 3 N 

HsM-^iO) 



r(li)=0 



$ 2 («) 



cr(u)=0 



a 1 



*i(u) 



o-(m)=0 



2 2 y^Jujo^u) cr|(u) 

a y^M ^>i(n)(g|(M) - <t?(m)) + 3$2(m)^(m)ct 2 (m)) 
2 i/o-^it) 



er(lt)=0 



(cr 3 (lt)cr 2i 3(M) - <7 2 (w)cr 3i3 (lt)) + <7 2 («)0 3 (u) - <7 3 (lt)(/> 2 (u) , 

</>(«) (cti(m)ct 3j3 (u) - ct 3 (u)cti :3 (m)) + ct 3 (m)0i(m) - ai(u)(j) 3 (u). 



Here <^(tt) and its derivatives were defined in (|59p . and the components u 2 ,w 3 of u are linear functions of 
ti , i 2 as described in (j!iO|) . These expressions give the solution of the Goryachev system for a general flow 
under the Hamiltonians H±, i/ 2 . 



5. Expansions of the functions x l ,y l and the Painleve analysis of the system. 

Apart from the formal cr-function solutions (|71l) it is important to know the complex singularities of 
the functions Ji(t),ji(t): their poles, order of branching, etc. Here we shall make connection with the 
Painleve analysis of |AvM989] . For simplicity we concentrate on the complex flow generated by the quadratic 
Hamiltonian Hi (time t = ti). As we have seen, this flow lies on the stratum C Jac(C), on which the 

coordinate u\ is a transcendental function of it 2 , u 3 and for which m 2 is constant along the flow. 

It follows from (fTTj) that most of the variables have poles when the w 3 -flow (i.e. the ii-flow) on 
crosses the substratum on which 0-3(11) = 0. Also, the right hand sides of (fTTj) may have branching 

only when the coordinate u\ (as a solution of a(u\, u 2 , it 3 ) = 0) ceases to be a locally meromorphic function 
of u 2 ,u 3 (and, therefore, of t\). As observed above, the implicit function theorem means this occurs along 
the 1-dimensional analytic subvariety Z = {<r(u) = 0, a\(u) = 0} C W 1 - 2 ^. (The solutions for J 2 ± iJ\ in 
(fTTj) have additional branching along due to presence of square roots.) 

To describe the local behavior of u\ — iti(u 2 ,u 3 ) in detail, choose a point Uq C Then cr(ito) = 

<t 3 (mo) = 0. Let Sui be the increments of the coordinates Ui such that Uq + Su C W^ 2 \ Then the following 
expansion holds 



(72) 



<t(u + Su) = cri(tto) 5ui + o- 2 («o) Su 2 + 



E 

l<ij<3 



= 0. 



For a generic Uq C and u 3 -flow (<5it 2 = 0) this implies 



(73) 



$ui = x(Su 3 ) 2 + 0((Su 3 ) 3 ) 



O-33{U ) 



2cri(u ) 

The above expansion does not hold for the points on W^ 1 ' — C with <Ji(u) = 0, i.e., at the points of ZnW^\ 

Proposition 5.1. The subvariety Z — {o~(u) = 0, o-\(u) = 0} has precisely 3 common points with W^ 1 : 
ZnWW = {r 1 ,r 2 ,r 3 }. 

Proof. In view of (|64j) . along the stratum one has o~\(u) = if and only if x\ = 0, which corresponds 
to the points r%, r 2 , r 3 . □ 
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Next, it follows from Proposition 14.21 that near u — ri, for the u^-Row (8u2 — 0), one has 



(74i a(r i + 8u) = a^(-- i (du 1 ) 



1 



p l 5 Ul 5u 3 + ^(Su 3 ) 2 ) + O (S Ul (6u 3 ) 2 ^ 



O((Su 3 ) 3 )=0, 

where, as above, p = exp(27n/3). Then the increment 8u\, as a function of du 3 , has 2 local branches 
(75) 8ux = \-p %i (-1 ± Vtj Su 3 + O ((Su 3 ) 2 ) . 

Indeed, substituting this into the expansion (|74|) . up to cubic terms, we get 

a(r t +Su) = ±^-(p^-l)(Su 3 ) 2 , 

which is zero for any i = 1,2,3. 

Note also that r x + r 2 + r 3 = in Jac(C) being the image of the divisor of x under the Abel map. 
Due to expansion ([24)) of the differentials flj near oo G C, a tangent vector to C at the origin 

u = 0, T W (1 \ is (0,0, 1) T . Hence the u 3 -flow is tangent to VF (1) at the origin. Next, in view of (|69j) . 
T r .W^ — (p l , 0, p 2l ) T , which means that the flow is also tangent to the projection of onto the (112, u 3 )- 
plane at r^. All these observations are depicted in Figure 5.1. 




FIGURE 1. A sketch of the stratum W^> C (solid line) and Z (dashed line) in the 

projection onto the (it2, U3)-plane, and the U3-flow on W^ 2 \ 



5.1. Expansions of the solutions near along the U3-flow. The order of poles of the variables Ji, ji, 
as functions of t\ (or u 3 ) depend on the nature of the intersection (for example, transversal or tangential) 
of the ii3-flow with C W^ 2 \ The expansions of the solutions ([TTj) in powers of u 3 may be found by 

using the corresponding expansions of the sigma-function and its derivatives near a point u$ G as 
well as the expansions (|73j) . ()75|) . It is more convenient however to find the expansions of the coordinates 
(xi,yi), (22,2/2) °f the points P\,P2 G S under the incomplete Abel map (1191) and then use the formulae 

GO). 

To do this, first note that, according to the definition of W* 1 ', when u belongs to \ {0}, one of the 
points Pi on C, say Pi, coincides with 00, whereas P2 remains finite. Now let T(t), \t — to\ < 1 be a complex 
analytic arc in such that T(to) H W' 1 ' — «o, and the projection of the arc onto C 2 — {u2 7 u 3 ) is a 



DYNAMICS ON STRATA OF TRIGONAL JACOBIANS 



17 



segment of a straight line U2 =const. Let D(t) = {Pi(t), P2(t)} be a divisor on C such that its Abel image 
gives r(t). Then P x (t ) = oc, P 2 (t ) = (x ,y ). 

Theorem 5.2. For a generic Uq = (wio, ""20, u 3o) £ W^ 1 ^, near i = t$ the coordinates of the points Pi, Pi 
admit the expansion 

Xl = UT. _ ^3 ~ — /,,„ _ „^2 + ° (("3 - "30) ) 

+ o((t~t )- 1 ) 



1 


3yo 


1 


("3 - W30) 3 


2rr (u 3 - 


- "30) 2 


27 1 


27 y 


1 


8\*-to) 3 


8 x (t 


-io) 2 


1 




1 


(«3 - "30) 4 


x (tx 3 - 


- U30) 3 


81 


27«y 


81 



(7(|) .'/i =- L - + — 7- — —— + ^11"; -";<>) Jl 



+ O((t-t )" 2 ) 



16(t-t ) 4 4a;o 16(t - i ) 3 



^2 = iC0 - 7^7— ("3 - U30) 2 + 0((u 3 - U 3 o) 3 ), 

2/2 = 2/0 + 0((u 3 - u 3Q ) 2 ). 
Next, near each point u$ = ri, i = 1, 2, 3, at which X20 — 0, i/iere are two expansions: 

1 1 3// 1 



X\ 



+ O((u 3 -u 30 )- 1 ), 



g 3 (u 3 -u 30 ) 3 2g 3 (u 3 - u 30 ) 2 

1 1 2p* 1 2 
2/1 = 774 7 71 T7 77 + °(( u 3 - "30j ), 



(77) ff 4 («3 - U30) 4 g 4 {u 3 - u 30 ) 3 

3 3 

X2 = -(-1 ± V3)p l (u 3 - u 3Q ) -g (u 3 - u 30 ) 2 + 0{{u 3 - -a 30 ) 3 ), 

2 p % 

3 

V2 = P 1 + jOe - M3P 4 ) (-1 ± \/3) 2 («3 - "so) 2 + 0((u 3 - u 30 ) 3 ), 

where g = (3 + y3)/2. 

The proof is given in Appendix. 

Theorem 15.21 implies that when the M3-flow crosses the substratum W^ 1 ) at a generic point tto, we have 

£22/1-^12/2 x , . 2/1-2/2 1 Lri , > 

= 1-0(1), = hO(u 3 -u 3 o), 

x 2 — xi 113- u 30 x 2 — xi 113 - "30 

(78) Xi+X2 = 7 rrr + 0((u 3 - u 30 )~ 2 ), x 1 X2 = - / — — -^ + O((u 3 -u 30 )~ 2 ). 

[u 3 - u 30 y (u 3 - u 3a y 

In the case of crossing at the points J7 g one has instead 

X2yi-X!y 2 . , y. i/i-Jfe 1 . nn s 
= O(l) + 0((U3 - U30)), = - + 0(1), 

X2-X1 X2-X1 g(u 3 ~u 30 ) 

! ' ! ' ; 1 ^/ 3p l (-l±73) 

.xi + x 2 = -^77 rr + 0(^3 - U30), xix 2 = 7-^7 77 + 0(1). 

g 6 (u 3 - u 3 or 2g J ( M3 _ Uw y 

We stress that in all the cases the above symmetric functions of Xi,yi, as functions of u 3 or t\, do not 
have poles with branching. However, they have finite branching along the intersection with Z C W^. 

Using (|71j) and the above expansions, one can estimate the leading terms of formal series solutions for the 
variables J\ ± zJ 2 , 72 ± «7i, J3, 73 near the poles. In the generic case one has 

73 = O (l/(f - to) 3/2 ) , J 3 = O (l/(t - t )) , 
(80) 72 + *7i=0(l), 72-*7i = 0(l/{t-t n ) 3 ) , 

Ji + 7 J 2 = O fl/(i - t ) 5/2 ) , Ji - zJ 2 = O f (i - io) 1/2 



18 H.W. BRADEN, V.Z. ENOLSKI, AND YU.N.FEDOROV 

and when the U3-flow crosses the substratum at r», 

73 = 0(1), J 3 = 0(l/(t -*«,)), 
(81) 72 + *7i =0(t- i ), 72 - %7i = O (l/(t - to)) , 

J x + 4 J 2 = o (l/(t - i„)) , J x - iJ 2 = 0(t - t ). 

These expansions correspond precisely with the leading behavior of formal series solutions of the Goryachev 
system obtained directly by Kovalevskaya-Painleve analysis. To observe this, first rewrite the Goryachev 
system (|3|) with the Hamiltonian Hi in the form 

d 4 
— (72 - 271) = -2i- J 3 (72 - 171) + 27 3 (Ji + 1J2), 

CLt\ o 

d 4 
-j— (72 + «7i) = J 3 (72 + 171) + 27 3 (Ji - 1J2), 
dti 3 



(82) 



d I T 1 M 2 T t 7 , 7 ^ 2 C7l +6 1/3 

— (Ji +iJ 2 ) = --JA-Ji +1J2) - g — (72 - m) + «a7 3 ' , 

— (Ji - iJ 2 ) = -J 3 H^i - W2) - g — 573— (72 +«7i) - Ja 7 3 



7 3 



d T 72 
- J 3 = -a- 



j| " 2/3 ' 

-77" 73 = + ^2X72 + «7i) - (^1 - «^2)(72 - «7l)- 

Then the corresponding formal Puiseaux (or Laurent) series solutions of (f82j) are of 2 kinds. According 
to terminology of |AvM989] . those depending on the maximal number of free parameters (here three, for 
example, the constants of motion h\, h 2 and a local coordinate on W^ 1 )) represent the principle balances of 
the solutions. Their leading behavior coincides with (|50)l . The series solutions depending on 2 or less free 
parameters (called secondary balances) correspond to (|81l) . 

Appendix. 

Proof of Proposition [4721 Let, as above, qi = (0, p l ) £ C, n = J J7. Since r, € W^ 1 ) and sc(&) = 0, 
Proposition 14. 1 1 implies that <7i(?7) = 03(77) = 0. Let £ be the local coordinate on C near 00 and x be such 
coordinate near qi = (0,//) G C. Introduce the functions 

U(0 = (Ui,U 2 ,U 3 )= [ O, and w(x) = {wi,W2,w 3 ) = [ fl, 

In view of the expansions (|24p . we obtain 

C/ ^-^-21^ 7 + --' 

(83) ^2 = -^ 2 -t^2£ 4 + t!m3£ 5 + --- , 

c/ 3 = -e + ^3e 4 + ^2C 5 + -- - , 

and, in view of (1691) . for i = 1 we have 

p 1 2 p 2 /i 6 - /i 2 3 1 4 1 (3p 2 /i 6 - 2/i 2 )M6 5 , rv 6% 
^=3^27 ~p X "iS^ 3 * +135 ~p X +0{X } ' 

P 2 1 2 p 2 M6 - M2 4 2 5 6 

w 2 = -x - — x - —pp 3 x +0(x ), 

6 36 p 45 



P 2 3 x 2 4, x " A^o 1 p6 ^2 p2 5 , 6\ 

™3 = ^- ^ - — p Me ^ - — M3 P a; + — 5 x° + 0{x°). 

3 27 36 135 p z 

The expansions for i = 2, 3 are obtained from the above by replacing p by p 2 and p 3 = 1 respectively. 
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Now note that for any i,(£Cwe have 7% + w(x) € W^- 1 ' and + w(x) + U(£) £ W^ 2 K Hence, in view 
of Proposition ^. 1[ 



E(x, := a{n + w(x) + [7(0 ) = 0, E 3 (x) = a 3 {n + w{x)) = . 
Substituting the above two expansions into the sigma-expansions (j65[) for i = 1, we obtain^ 



. , / \ 1 4 1 2 1 

£(a;,f) = I g ct 2 p+ — cr 3 , 3 p + Y8 CTl - 1 ^ + 9 fJl ^ 



-g(o-3,3P 2 +ffl,3P)^ + (~^ CT2 + ^ CT 3,3)f 2 

1 2 1 3 1 1 4 \ , 2 

^ffl,l,3P +g01 )3 , 3 p + -C 2 ,3/0+ ^°"3,3,3P I 4 ^ 

1 2 1 1 1 2 \ /l 1 \ - , 

-gCT 2 ,3/0 + g CT 1:3j3 p- -cr 12( o+ - (73,3,3 P I ? X + I - CT 2 , 3 - - <7 3 , 3 , 3 I 4 + 0(£ x) 



and 



/ n 1 2 \ ^ 1 4 1 2 1 3 1 \ 2 n , 3 , 

E 3 (x) = I g CT1, 3 p + - CT 3j 3 p I X + I — CT 3 , 3, 3 p + — (T1,1,3P + g 171, 3 , 3 p + - <7 2 , 3 P I X + 0(0! ). 

Then, upon equating the coefficients of the expansions E(x, £), S 3 (a;) to zero, we get a system of equations 
for the coefficients it- k . (Note that the indicated coefficients of E 3 (x) are also coefficients of E(x, £), so 
they, in fact, do no bring new conditions.) Solving it, we obtain the first group of relations (l66l) . 
To find the relations (|67|) . we proceed in the same manner as in ([53]) and consider the limit 

,. oifo +u) + pV;(7-; + u) 
Q% = hm 



u^o a^n + u) 

Applying l'Hopital's rule, we get 

q _ gi(r» + u 3 e 3 ) + p* gj(r^ + tt 3 e 3 ) _ ^ Ji 3 (r^ + M 3 e 3 ) + p' cr 33 (r^ + M 3 e 3 ) 
ii 3 ->-o ^(ri + uses) U3-S-0 cr 3 (ri + w 3 e 3 ) 

= Cl33(^) +p' 0-333 (^) 
0-33 

On the other hand, 



e 3 = (0,0,l) J 



,. 01 0*t + U2e 2 ) + P cr 3 (r.i + U 2 e 2 ) O^fa) + p 1 <7 2 z{ri) /n 1 nVT 

hm ■ r = — — , e 2 = (0,1,0) 

«2^o cr(r J + u 2 e 2 ) cr 2 (r;) 

CTi(r l + mei) + p l <7 3 (ri + mei) 011 fa + inei) + p* oi 3 fa + uiei) 

lim ; - = lim 



"i->o cfa + uiei) ui-Hi (Ti(ri + uiei) 

^lnfa) +p* o- n3 fa) T 

o-nfa) 

Comparing the above limits and using (|66|) . we get 

0-133(^3) + 3pV 23 (r 3 ) = CT12 (r 3 ) + pV 23 (r 3 ), 

_ j Olllfa) + p' 0-113(^3) _ ^12 fa) + P' 023 fa) 

P 2cr 2 (r 3 ) o- 2 (r 3 ) 

These relations are compatible with a vanishing of the leading coefficients of the expansions E(x,£),E 3 (x) 
if and only if relations (|67l) hold. This proves the proposition. 



Here we omit the index (i) in the cr-derivatives. 
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Proof of Theorem \5.2[ It is sufficient to study the expansion of the map (fT9"f near the point {oo, (xo,yo)} £ 
C x C. Again let £ = a; -1 / 3 be a local coordinate of Pi S C near oo and n = x — x be such a coordinate of 
P2 near (a:o,yo)j ^o 7^ 0- The expansions of the differentials near (a;o,yo) are 

a; 



2 = 
^3 



dG/dy(x,y) 

y 



dx = (<p 2 o + + 0{k 2 )) dn, 
dx = (1^30 + ip 3 iK + 0(k 2 )) ck, 



dG/dy(x,y) 

where ip 2 o, ■ ■ ■ ,<f3i are certain nonzero constants depending on xq and the coefficients of the curve. Note 
that 4^ = In view of this and of the expansions of Qj in (|24|) . the differential of the map (|19|) reads 



(84) 



£+^£ 3 + 0(£ 4 ) ^20 + W + O( K 2 )' 



f 



dw 2 
dw 3 



Taking into account that <p20,<P30 ar e non-zero, we may invert the above matrix expansion yielding (up to 
quadratic terms) 



(85) 



'30 
'20 



'30 



- <E>K 

^21 



f 



>!0 



V02O 0l O ( '-]u "20 



du 2 
du 3 



$ = 



31 920 



»21<P30 



>20 



|zcft and du 2 = 0. Then (|55|) gives a system of 2 ODEs with respect to 



For the ^-flow we have du 3 
u 3 ~ u 30 ( or t ~ to)- Applying the condition £(0) = 0, k(0) = 0, we find initial terms of the series solutions 



£ = u 3 



"30 + ~{u 3 ~ u 30 ) 2 + 0((tt 3 - "so) 3 ) = §*(t - t ) - - ^o) 2 + ■ 



f 



202 



■(«3- "30) 2 + 0((lt 3 - U 30 ) 3 ) 



1 2 s 
2020 3 ; 



(t - *o) 2 + 



Substituting this into the expansion (|22[) and taking into account f|54|) . we get the expansions (|76|) . 

The above argument fails to work only when i^^o) = ?i = (0, p 1 ) and 02o(<Zi) = 0, i.e., when the differential 
relation (|84|) cannot be locally inverted. In this case we will use the expansions of the differentials fii near 
qi given by (|69l) and replace (|84|) by the differential of the Abel map (|20|) taking the part with Oi , O3 and 
U\, 113. Let i be a local coordinate on C near = (0, p 1 ). For a — 1 the expansion of the differential is 



(86) 



e + + o(e) f - ^pVe - ^ + o^) 



i 
3 



PeP 



x 2 + 0(x 3 ) 



3 ' ' ~ v ' ' 3 9 

The latter is invertible and, up to cubic terms in £, x, gives the expansions 

'-P - P 2 £ + 3^ (-Me + PP2) z 2 - £ 2 1 + £ p + p 2 £ 2N 



(87) 



+ 3£+(2p 6 p 2 - m )x 2 + 3p£ 2 



3£ 



-3£ 2 



\du 3 



According to (f75|) . for i = 1 one has diii = vdu 3 , v — |(— 1 ± \/3)p 2 - Then the above expansions lead to the 
following system of 2 ODEs with respect to u 3 — u 3 q 



^ = g{ i +Pi+P 2 e 

du 3 

£L = ^ _ ^£ _ 3ff £ 2 + (2p 2 p 6 - p 2 ),* 2 + 0(£* 2 ), 
«u 3 p p 



T ^(PM2-P 6 )a; 2 + C>(£ 3 ) 
3p 2 



where 



= 1 -/»«/ = l-i(-l±>/3) = 



3T V3 
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Applying again the conditions £(0) = 0,x(0) = 0, we find the series solutions with initial terms 

6 = g(u3 - U 3Q ) + ^7f("3 - U 30 ) 2 + O ((u 3 - U 30 f) , 

x = — (u s - u 30 ) - —{u 3 - u 30 ) 2 + O ((u 3 - u 30 ) 3 ) , 
P P 

which, in view of the expansions (|22|) . (|68|) . give (|77|) . The expansions for i = 2,3 are obtained by replacing 
p by p 2 and by p 3 = 1 respectively. □ 
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